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Abstract 

Nonlinear propagation of electrostatic modes in ultrarelativistic dense elelectron-positron gravito- 
plasma at the polar cap region of pulsar magnetosphere is considered. A nonlinear Schrodinger 
equation is obtained from the reductive perturbation method which predicts the existence of 
Langmuir dark solitons. Relevance of the propagating dark solitons to the pulsar radio emission is 
discussed. 
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I. Introduction 

Pulsars are celestial sources that believe to be rotating 
neutron stars producing light-house like beams of radio 
emissions from the magnetic poles. As shown by Goldre- 
ich and Julian (1969) the rotating magnetic dipole pro- 
duces a quadrupole electric field whose component paral- 
lel to the open magnetic field lines at the poles extracts 
particles very effectively from neutron star surface and 
accelerates them to highly relativistic energies. Thus, 
the magnetosphere is filled with plasma which shields 
the electric field. Complete shielding is established when 
the net charge reaches ncj - the Goldreich - Julian charge 
density. The Lorentz factors of the accelerated particles 
reach about 10^ and they emit hard curvature radiations 
that propagate at a sufficient angle to the magnetic field, 
so that significant pair production of electron-positron 
can occur (Erber 1966). It is commonly accepted that the 
newly created particles produce more pairs by emitting 
energetic synchrotron or curvature radiation. As a result 
an avalanche of secondary particles populates the magne- 
tosphere with densities W^noj (Ruderman and Suther- 
land, 1975) . Here, we extend our earlier research on pul- 
sar microstructure, soliton formation, wakefield accelera- 
tions, gravitational waves , and growing modes (Mofiz, et 
al. 1985-2011) to account the pair ultrarelativistic pres- 
sure and the gravity 

The paper is organized as follows. Section II describes 
the fluid model of the dense ultrarelativistic electron- 
positron plasma under gravity. Considering a Lorentz 
invariant frame moving with group velocity of the wave, 
a Nonlinear Schrodinger Equation (NLSE) is derived us- 
ing the reductive perturbation method (Gardner and 



Morikawa, 1960). A linear dispersion relation is obtained 
showing the existence of Langmuir waves under gravity 
and with ultrarelativistic temperature for wave propa- 
gation. The solution of NLSE shows the generation of 
Langmuir dark solitons. Results are discussed in Sec. 
III. Finally, Sec. IV concludes the paper. 



II. The Mathematical Model 

We consider two-fluid magnetohydrodynamic (MHD) 
equations to describe the electron-positron plasma in the 
pulsar magnetosphere of the neutron star. The equa- 
tions are the usual continuity and momentum balance 
equations for the plasma species, electrons and positrons, 
supplemented by the Poisson's equation for electrostatic 
wave propagation . Thus, the required set of equations 
are as follows (Mofiz and Ahmedov, 2000): 
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where, 7s = 1/a/1 — v'^/c^, Hs, Vg, and Qs are respec- 
tively the particle number density, particle velocity, and 
particle charge of plasma species s; qs ^ —e ioi s ~ e 
(electron) and = -\-e for s = (positron); cp is 
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the electrostatic potential, Tg is the Schwarzschild ra- 
dius of the neutron star; m is the electron/positron 
mass; e is the absolute value of the electronic charge. 
In Eq. (2), the pressure Ps is given by the expression 
for the ultrarelativistic pressure (Chandrasekhar, 1938): 
Pa = noksT {us/nof'^^ , where = Tp = T has been 
assumed; fc^ is the Boltzmann constant, no is the equi- 
librium particle number density. 

In the polar cap region of the pulsar, we consider 
^ = 0, V = zd/dz, Vs = VszZ and adopt the follow- 
ing normalization of different quantities: z — )• zujpe/c, 



t tujpe, ris njno, Us v^Jc, <j) -4 e<j)/mc'^, 
gt — > kBT/mc^, Vg — > rgUjpe/c, and Vg — > Vg/c, where 

1/2 

* is the electron plasma frequency. 



c is the speed of light. To study the nonlinear dynam- 
ics, we consider the following stretched coordinates in the 
moving frame (Melikidze et al. 2000): 



and 
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where, Vg is the group velocity, 70 = 7,. = 7p is the aver- 
age relativistic Lorentz factor and is given by the follow- 

ing expression: 70 = (l — t'g) , e is a small quantity, 
the perturbation parameter with e < 1. With the trans- 
formations, given by Eqs. (4) and (5), the derivatives 
d/dt and d/dz transform to d/dt — > d/dt — €'joVgd/d£^ + 
e^Jod/dr and d/dz — > d/dz + ej^d/dS, — e^^QVgd/dr 
respectively. With the above considerations, the set of 
equations, Eqs. (l)-(3) take the following forms: 
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where, the factor g{z) = (1 — Tg/z) ' accounts for the 
gravitational effect. Now we expand the quantities rig. 
Us, (f) as 



= 1 + e^Hso 
00 



1=1 



Us = e^uso + ^ e' (usie'^^'''-^'^ + ^^e-^'C^^""^*)) ,(10) 
1=1 

00 
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Here, 

(n,o, Wso, 00 , nsi, Usi, = tK^) + eA(2) + e^^^^) ^ 

are functions of stretched coordinates (^,t). 
II. A. Linear Dispersion Relation for the Lang- 
muir Wave 

Now considering \\/ g{z) ■ dg{z) / dz\ « k for the first 
harmonic {I = l)in the first order (e = 1 ), we have the 
following equations for the first-order quantities: 
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Eliminating u^^^ from Eqs. (12) and (18), we obtain the 
following equation relating n^^^ and (^^^^ : 
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from which we obtain 
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Using Eq. (16) into Eq. (14), we obtain the following 
linear dispersion relation: 



which in the dimensional form is 
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The group velocity Vg = dui/dk is 
obtained from the linear dispersion relation Eq. (17) as: 
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The same expression for Vg is also obtained from the 
compatibility condition and shown in the Appendix A. 



The group dispersion is found to be 
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Eq.(17) represents the dispersion relation for Langmuir 
waves in ultrarelativistic e, e+ plasma under gravity. Pair 
production in the polar cap region of pulsar magneto- 
sphere occurs through curvature radiation which happens 
for £\\ » rrieC^, where £|| is the energy of electron along 
the magnetic field. It is estimated that for cascade gen- 
eration of pair plasma 70 ~ 10*3 - lO^fy ~ lO^^ - lO^^eV 
(Beskin et al. 1993). Here, 70 = f|| /mgC^ ,then consid- 
ering 5[| = ksT, we find 70 = ksT/meC^ = gt- Using 
Eqs.(17),(19) and Eq.(20), we perform an analysis of the 
dispersion relation, group velocity and group dispersion 
of Langmuir waves at the ultrarelativistic temperature 
of the e, e"*" plasma under gravity. The analysis is shown 
graphically in Fig. 1-4, respectively. 

II.B Nonlinear Evolution Equation for the Lang- 
muir Wave 

Finding the zeroth harmonic and second harmonic of the 
second order quantities in terms of the first harmonic of 
the first order quantities and using these into the first 

harmonic of the third order quantities , we easily obtain 
the following NLSE for the evolution of the potential 



.da „9^a ^, ,2 



0, 



(21) 



where a = 



and the coeflScients P and Q are given 



by the following expressions: 
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FIG. 1: The variation of the normalized frequency w of the 
linear Langmuir wave with respect to the normaUzed pump 

wavenumber k for dilferent values of plasma parameters : = 
1, = 2, 70 = ut = 10^ with the corresponding temperature 
T = 5 X 10^^ K . 
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FIG. 2: The variation of the normalized group velocity Vg 
of the linear Langmuir wave with respect to the normalized 
pump wavenumber k for different values of plasma parameters 
: rj = 1, z = 2, 70 = (Jt = 10^ with the corresponding 
temperature T = 5 x lO^'^ K . 
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FIG. 3: The variation of the normalized group dispersion v'g 
of the linear Langmuir wave with respect to the normalized 

pump wavcnumber k for different values of plasma parameters 
: Tg = 1, z = 2, 7o = (Tt = 10*' with the corresponding 
temperature T = 5 x 10^^ K . 
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FIG. 4: The variation of the normalized frequency w of the 
linear Langmuir wave with respect to the normalized distance 
z for parameters : Vg = 1, k = 1000,7o = ar = 10® with the 
corresponding temperature T = 5 x 10^^ K . 
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Here, the coefficient P can be written as 
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represents the effects of ultrarelativistic temperature 
and gravity , neglecting of which we recover the usual 
results. The coefficients P and Q appearing in the NLSE, 
given by Eqs. (22), (23)are known as the dispersion and 
nonlinear coefficients, respectively. The signs of P and 
Q determine whether the slowly varying wave pulse is 
stable or not (Lighthill condition; Lighthill, 1967). If the 
signs of P and Q are such that PQ < 0, the wave pulse is 
modnlationally stable and the corresponding solution of 
the NLSE is called the dark soliton. On the other hand, if 
PQ > 0, then the pulse may be modnlationally unstable 
and the solution of the NLSE in this case is called the 
bright soliton. Graphically, we study the nature of P 
and Q for continuous values of the wave number k with 
particular values of plasma parameters, which are shown 
graphically in Fig.5-6, respectively. 
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FIG. 5: The variation of the dispersion coefficient P of the 
linear Langmuir wave with respect to the normalized pump 
wavenumber k for different values of plasma parameters : Vg = 



1, 2; = 2, 70 = (7T = 10" 

T = 5 X lO^'^ K . 



with the corresponding temperature 



From the graphical analysis, we find that PQ < 0. 
Then applying the standard technique, and by taking 
C'(C) ■'") = ~ ^t)], the following solution of 



5 



-1x10'^ 



O'-2xl0" 



-3x10" 



-4x10 



13 



500 



1000 1500 2000 



2.X10- 



..5x10- 



1.x 10" 



5.x 10-1" - 




FIG. 6: The variation of the nonlineax coefficient Q of the 
linear Langmuir wave with respect to the normalized pump 

wavenumber k for different values of plasma parameters : Tg — 
1, 2; = 2, 7o = aT = 10^ with the corresponding temperature 
T = 5 X 10^5 K . 



FIG. 7: Dark soliton in ultrarelativistic e, e"*" plasma at the 
polar cap region of pulsar magnetosphere. The parameters 



arc : rg = 1, z = 2, fc = 1000,0 = 0,K = 1„ 70 = ctt 
with the corresponding temperature T = 5 x 10^^ K . 
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the NLSE, Eq. (21) (Mofiz,2007)is easily obtained: 

T 1/2 



a(^, t) — Go tanh 
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}iere,K£^ — fir is the modulation phase with /!'(<< A;) 
and 0(<< uj), respectively. Eq.(30) represents a dark 
soliton with amplitude 
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and width 6 
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, respectively. 



The dark soliton (Eq.(30)) in the ultrarelativistic e,e-\- 
plasma is shown graphically in Fig. 7. 



III. Results and Discussion 

In this section, we analyze the linear dispersion as well 
as the nonlinear Langmuir dark soliton in the pulsar mag- 
netosphere. Eq.(17) shows that the Langmuir frequency 
depends on ultrarelativistic temperature and it is red- 
shifted due to gravity near the Schwarzchild radius. Sim- 
ilarly, the group velocity and group dispersion, shown by 
Eq.(19) and Eq. (20), arc also depend on temperature 
and gravity. 

For numerical appreciation of the dark soliton ,we 
consider the two cases of ultrarelativistic temperatures: 
6 X IQ^^K - 1.8 X 10"^ (Crab pulsar) with the corre- 
sponding energies 5 — 15MeV (Nanobashvilli,2004) and 
5 X 10^^.^ — 5 X 10^^ K (x-ray pulsar) with the correspond- 
ing energies lO^^ _ iQ^s^y (Bcskin et al., 1993). 

The solution of the NLSE (Eq.(21)) is a stable dark 
soliton (Eq.(30))whose amplitude and width depend 
on temperature. The amplitude is increased and the 
width is decreased with the increase of ultrarelativistic 



temperature. Thus, stable spiky Langmuir solitons are 
possible in the ultrarelativistic electron-positron plasma. 



IV. Conclusion 



To summarize, we have investigated the nonlinear 
propagation of electrostatic modes in a dense ultrarela- 
tivistic electron-positron gravito-plasma at the polar cap 
region of pulsar magnetosphere. A multiscalc perturba- 
tion analysis of the fluid equations shows that stable dark 
Langmuir solitons are produced due to the balnce of dis- 
persion and nonlinearity in the wave propagation. As the 
amplitude of the soliton increaes and width of the soliton 
decreases with the increase of ultrarelativistic tempera- 
ture, so spiky stable dark Langmuir solitons may prop- 
agate along the open field lines of the pulsar magneto- 
sphere, which may have some relation with pulsar radio 
emission and its microstructure. 
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Appendix A: The compatibility condition 

It can be shown that the 1st harmonic of the 2nd-order 
electron and positron densities can be found to be 
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After finding n^j^ — rip^ and substituting it in the fol- 
lowing Ist-harmonic of the 2nd-order part of the Pois- 
son's equation: 



gizf' 70 ^n,, n^,)+ ^^^^ 



0, 



we obtain the following compatibility condition: 
_ 4aTg(z) k 

which is exactly the same as the expression of the group 
velocity, Eq. (19), obtained by differentiating oj with 
respect to k from the linear dispersion relation, Eq. (17). 
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